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ABSTRACT 
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An experimental apparatus was designed, constructed and 
instrumented in an effort to systematically and carefully 
study the condensation heat- transfer coefficient on a single, 
horizontal tube. A smooth, thick-walled copper tube of 
length 133.5 mm, with an outside diameter of 15.9 mm and 
an inside diameter of 12.7 mm was instrumented with six 
wall thermocouples. The temperature rise across the test 
section was measured accurately using quartz crystal 
thermometers. The inside heat- transfer coefficient was 
determined using the Sieder-Tate correlation with leading 
coefficient of 0 . 029 . Initial steam side data were taken 
at atmospheric pressure to test the data acquisition/ 
reduction computer programs. 
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I. INTRODUCTION 



A. SIGNIFICANCE 

Experimental studies have shown that enhancement techniques 
applied on the steam side of horizontal condenser tubes may 
significantly reduce the thermal resistance and thus reduce 
the overall size and weight of the condenser. Several work- 
shops and conferences on heat exchangers have addressed 
specific areas in need of research {Refs. 1,2,3). Recently, 
there have been many investigations into possible alterna- 
tives for enhanced heat-transf er surfaces in the design of 
heat exchangers {Refs. 2, 4, 5, 6, 7}. Webb {Ref. 8} in 1980 
published an excellent review of enhancement methods as 
applied to condensers. 

Enhancement techniques can generally be divided into 
two major categories: water-side enhancement and steam-side 

enhancement. Water-side enhancement is covered extensively 
by Bergles & Jensen {Ref. 9). Various steam-side enhance- 
ments include low integral fins, surface coatings, fluted 
tubes and roped tubes. 

The goal of any enhancement technique is to lower the 
thermal resistance to the point where the technique becomes 
economically feasible. Althouth Bergles & Jensen {Ref. 9) 
show that extensive reduction in thermal resistance can be 
obtained on the water-side, it is felt that a significant 
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portion of the overall thermal resistance can be reduced by 
enhancement techniques on the steam-side. Both methods must 
be examined independently because the best enhancement for 
one side may not be the best for both sides. 

Several research efforts have been undertaken at the 
Naval Postgraduate School to model steam surface condensers 
and analyze their thermal performance numerically {Refs. 
10,11,12}. Based on numerous computer runs used to analyze 
different combinations of both inside and outside enhancements, 
it is felt that further experimental research should be 
conductdd on steam-side enhancement techniques. 

B. OBJECTIVE 

The objective of this research effort was to obtain 
smooth-tube baseline data necessary to act as the standard 
for future enhancement techniques. Both the inside and 
outside heat-transfer coefficients had to be identified 
and all experimental methods used to obtain and reduce the 
recorded data had to be well documented so that future 
testing using enhanced tubes would be performed under identical 
conditions. 
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II. DESCRIPTION OF EXPERIMENTAL APPARATUS 



A. OVERVIEW OF SYSTEM 

The experimental apparatus used in this investigation, as 
described in xReference 13, was constructed of stainless steel 
and 152.4 mm (6 in) diameter Corning Pyrex glass piping. An 
overall schematic of the system is shown in Figure 2.1 and 
a photograph is shown in Figure 2.2. A boiler section 
304.8 mm (12 in) in diameter and 487.7 mm (19.2 in) in 
height, with a 304.8 mm (12 in) to 152.4 mm (6 in) reducing 
section, provides steam for the system. Ten Watlow 
4000-Watt, 480-Volt stainless steel immersion heaters are 
used to provide input power. Each heater is 15.9 mm (5/8 in) 
diameter and 279.4 mm (11 in) long. All heaters are 
connected in parallel. 

Stean travels up the glass piping to a height of 3.1 m 
(10 ft) and makes a 180-degree bend followed by a 1.52 m 
(5 ft) straightening section prior to passing through the 
condenser test section. Two helically-wound, water cooled 
coils of 9.5 mm (3/8 in) diameter copper tubing act as a 
dumo condenser to remove all excess steam. A 25.4 mm (1 in) 
diameter stainless steel tube returns condensate by gravity 



to the boiler. 



B. CONDENSER TEST SECTION 



The condenser test section is made of stainless steel, 

152.4 mm (6 in) in diameter and 457.2 mm (18 in) in length. 
Figure 2.3 shows a schematic and Figure 2.4 shows a photo- 
graph of the test section. Steam is condensed on a single 
horizontal tube. Two types of smooth test tubes were 
constructed of high-grade copper to provide baseline data for 
inside and outside heat- transfer coefficients. An instrumented 
tube, constructed as shown in Figure 2.5, and an uninstrumented 
tube of the same dimensions were used to obtain initial data. 
The unins truraented tube was constructed identical to the 
instrumented tube but there were no channels machined for 
thermocouple installation. The active condensing length 
was 133.3 mm (5.2 5 in) with a suitable correction applied 
for the "fin” effect created by both ends of the tube not 
being in direct contact with the steam on the outside, 
but cooled by water on the inside. The fin-effect calculations 
are described in Appendix A. 

The test tube is held in place by nylon holding rings as 
described in Reference 13. An additional 0-ring was placed 
inside the inlet section to prevent cooling water from 
leaking into the grooves of the instrumented test tube. 

Nylon was used to minimize axial conduction from the test 
section. A nylon mixing chamber described in Reference 13 
was placed in the outlet side prior to measuring the temper- 
ature to ensure thorough mixing of the coolant. A 114-3 mm 
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(4.5 in) diameter double glass window, heated by hot air 
to prevent fogging, was used to provide visual observation 
of the condensing process. 

A 12.7 mm (0.5 in) port at the dump condenser section 
was used as an overpressure relief line and boiler fill 
line. An in-line 0.108 MPa (1 psig) relief valve was used 
to allow gases to escape while warming up the system. 

A 12.7 mm (0.5 in) vacuum valve was placed in a branch line 
off of the relief line to allow filling of the boiler as 
described in Appendix B. All connections were assembled 
with Swagelok fittings sealed with Teflon ferrules to ensure 
leak tighness. 

1 . Cooling Water System 

A closed loop cooling water system was designed and 
constructed to circulate cooling water through the test 
section. A 0.001 cu-m (2.8 cuft) stainless-steel tank was 
used as a reservoir which fed the suction end of a close 
coupled centrifugal pump rated at 0.0022 m / s (11 GPM) at 
12.19 m (40 ft) of head. A 19.1 mm (3/4 in) diameter bronze 
needle valve was used to regulate the flow which was 
monitored by a rotameter with a full-scale range of 6.94 
E-4 m^/s (11GPM). See Figure 2.6 for details. The maximum 
flow rate through the test tube was 4-57 m/s (15 fps) at 
an inlet temperature of 15.6°C (60°F). The centrifugal 
pump needle valve arrangement provided a flow rate accuracy 
of 0.5 percent on a scale with 100 divisions. After leaving 
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the condenser test section, the water was cooled back down 
to the 15.6 C ( 60 F ) bath temperature by means of a 11.73 kW 
( 40 , 000-BTU/Hr ) direct expansion heat exchanger. The heat- 
exchanger refrigerant dissipates the heat to the atmosphere 
through a 10.55 kW (3 ton) air-conditioning unit. The 
chiller section is illustrated in Figure 2.7. 

The closed loop cooling water system was chosen in 
order to more closely monitor the water flow rate and the 
inlet temperature to the test condenser tube. The reservoir 
is easily filled by use of the installed fill line and 
drained by use of a siphon tube. 

2. Test Condenser Tubes 

Initial smooth tube data were taken using an instru- 
mented copper tube as illustrated in Figure 2.8. Construction 
of the instrumented tube was carefully performed in accordance 
with Figure 2.5. The instrumented tube must be installed 
according to the instructions as outlined in Appendix D to 
prevent damage to thermocouple wires and leakage from the 
system. Each wall thermocouple should be checked for proper 
operation prior to taking any data. 

C. DUMP CONDENSER SYSTEM 

A dump condenser is installed to condense excess steam 
from the test section. Cooling water is supplied from tap 
water through two parallel, helically-wound, 9.5 mm (3/8 in) 
diameter, copper coils. One smaller diameter coil is placed 



inside the other, while both are 457.2 mm (18 in) high. The 
tap water flow rate is controlled by a 19.05 mm (3/4 in) 
diameter bronze needle valve and is monitored using a 
rotameter. System internal pressure can be controlled by 
reducing the flow of cooling water to the dump condenser. 

For a given boiler power input, a dump condenser rotameter 
reading will correspond to a system pressure. 

D. VACUUM SYSTEM 

Two operational vacuum systems are available to maintain 
sub-atmospheric conditions within the apparatus; an air 
ejector, powered by 1.1 MPa (160 psig) air supplied by the 
Mechanical Engineering Department reciprocating air corn- 
pressure, and a portable mechanical vacuum pump (Figure 2.6). 
Either system may be used to initially draw down the pressure 
but only the air ejector may be used while the system in in 
operation due to the lack of a nitrogen cold trap on the 
mechanical vacuum pump. The vacuum system inlet valve can 
be used to throttle the air ejector to maintain any system 
pressure above 13.8 kPa (2 psia). This system is also 
useful in removing noncondensable gases from the apparatus. 
The sequential drawing down of vacuum on the system while 
steaming, shutting the vacuum inlet valve and raising the 
pressure to about 13.8 kPa (2 psia) by adjusting the dump 
condenser throttle valve was found to be an effective 
method of removing noncondensable gases. This procedure 
should be repeated about three times. 
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The system will not maintain vacuum due to cracked silver 
soldered connections at the base of the dump condenser 
section. The copper tubing used as a dump condenser pene- 
trates the stainless steel base and was sealed with silver 
solder which cracked, producing small leaks. 

E. INSTRUMENTATION 

1 . Heater Power Control 

Power input to the apparatus is regulated by a 
voltage sensing circuit. Line voltage (440-VAC) is fed into 
a differential input precession voltage attenuator which 
divides the line voltage by one-hundred. This stepped- 
down voltage passes through a True Root Mean Square converter 
stage on which the integrated period has been reduced to 
about 1 ms. The output of the TRMS converter is buffed and 
compared to a reference voltage provided by a potentiometer 
mounted on the front panel. See Figure 2.9 for a general 
view of the control panel. The comparator output is fed 
to the control input of a Halmar silicon-controlled rectifier 
power supply which provides the voltage applied to the heater 
elements. The buffered output of the TRMS converter is 
amplified and filtered to provide a normalized voltage 
output proportional to the actual voltage output of the 
power supply which is monitored by the data acquisition/ 
reduction system. A line diagram is provided in Figure 
2 . 10 . 
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2. Pressure Measurement 



A U-Tube mercury manometer, calibrated in millimeters, 
is installed to measure system pressure. The pressure tap 
is located at the test section and is tilted upward to avoid 
steam from collecting in the manometer (Figure 2.4). A 
vacuum valve is installed in the pressure line to prevent 
condensation from accumulating in the manometer. A small 
amount of water eventually collects in the manometer, but 
this discrepancy is taken into account by the data acquisition/ 
reduction software with manual input to the keyboard of 
mercury height followed by water height. 

3 . Temperature Measurement 

Various temperatures are monitored throughout the 
system to include: 1. test section cooling water inlet 

and outlet, 2. steam, 3. condensate return, 4. ambient, 

5. test tube wall and 6. cooling water reservoir. The 
temperature rise across the test section is measured by a 
Hewlett-Packard quartz crystal thermometer accurate to 
0.04°C. The location of the temperature probes are shown in 
Figure 2.4* The cooling water reservoir temperature is 
monitored by a direct readout alcohol-in-glass thermometer. 

All other temperature measurement is accomplished by 
copper-constantan thermocouples: two for the steam, six for 

the test tube wall and one for the condensate return. Each 
temperature measurement, except for the cooling water 
reservoir, is read directly by a Hewlett-Packard 3497A 
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data acquisition system, which is controlled by a Hewlett- 
Packard 9826 computer. 

During operation, the six wall thermocouples 
experienced temperature fluctuations on the order of - 10 uV. 
Therefore, each thermocouple was scanned for ten seconds 
and the ten readings were averaged to obtain a more accurate 
measurement. All other thermocouples maintained a steady 
readout with only 1.0 uV fluctuation. The quartz crystal 
thermometer appear ed to operate very well with no fluctuations. 

Two of the thermocouples were calibrated. One 
thermocouple was made from the beginning of a spool of 
copper-constantan wire and one was made from the end. 

All other thermocouples were fabricated from wire on this 
spool. It was assumed that the properties of the copper 
and constantan do not change along a given section of wire 
for any given spool. Both thermocouples were calibrated 
by the method described in Appendix E. 
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Figure 2.1 Schematic of Test Apparatus 




Figure 2.2 Photograph of Overall System Showing 
Data Acquisition System 
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Figure 2.3 Schematic of Test Section 
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Figure 2.4 Photograph of Test Section 




Q : 
<n 
m 





27 




! 




Control Valve 



gp® M£ML 






mKHHb- 

I Cooling 
I Water 
1 pumD 






i;/>: 



Cooling Water 
Reservoir 



Ro tamet er 



Figure 2.6 Photograph of Cooling 
Water System 
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Figure 2.7 Photograph of Overall System Showing 
Boiler and Air Conditioning Unit 
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Figure 2.8 Photograph of Instrumented Test Tube 
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Figure 2.9 Photograph of Control Panel 
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Figure 2.10 Power Controller Circuitry 



III. SYSTEM OPERATION 



After a thorough cleaning using the procedures described 
in Appendix F, a test tube may be installed by the steps 
outlined in Appendix D. Each system start-up should be 
accomplished by the procedures described in Appendix B. 

Once the system has reached steady-state conditions (about 
30 minutes), data taking for a given test tube may commence. 
During each data run, a frequent check of the test tube 
should be made to ensure the desired condensation mode: 
either dropwise or filmwise. When the tube cleaning procedures 
of Appendix F are utilized, filmwise condensation conditions 
will last for about four hours, followed by mixed filmwise 
and dropwise conditions. Each data set consisted of eight 
cooling water velocities ranging from 1.37 m/s to 3.84 m/s 
through the test section. At a given velocity setting, the 
flow was allowed to steady out for two minutes before 
recording temperatures. The computer signalled by both 
audio and visual means when the system was ready for a new 
data point. An entire data set was recorded in 40 minutes. 
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IV. DATA ACQUIS IT ION /REDUCTION 



A. DATA ACQUISITION AND STORAGE 

A Hewlett-Packard (HP) 3054A automatic data acquisition/ 
control system was used to monitor temperatures from the 
thermocouples and the quartz crystal thermometers. Figure 4.1 
shows a photograph of the system. This system included an 
HP 3497A data acquisi tion/control scanner and an HP 3456A 
digital voltmeter with a resolution of 1 00-nano vol t s . An 
HP 9826A computer, interfaced with the data acquisition/ 
control system, was used as a controlling unit and provided 
storage for data. Information was entered through the key- 
board to prompt the data acquisition/control unit to auto- 
matically scan each channel. Channel assignments are listed 
in Table 1 . These raw data were then transferred to a 
computer disk under a user-specified file name for later 
reduction to usable form. The ability to store raw data 
directly enabled these data to be reduced at any time and 
allowed flexibility for changes to da ta-reduction software. 

3. DATA REDUCTION 

Following data acquisition for each data point, results 
were computed according to the stepwise procedure outlined 
in the next section, and were then printed on an HP 6771 G 
thermal printer. Appendix G shows an example of a represen- 
tative data run. Results were also stored in a user-specified 
plot data file for subsequent plotting using the program PLOT. 
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C. STEPWISE SOLUTION PROCEDURE 
1 . Program SIEDER 

a. Compute average cooling water temperature. 

b. Compute average wall temperature. 

c. Compute cooling water velocity. 

d. Compute mass flow rate of cooling water. 

e. Compute heat transferred to cooling water. 

f. Compute average inside wall temperature. 

g. Compute LMTD. 

h. Assume fin efficiency. 

i. Compute inside heat- transfer coefficient. 

j. Compute Nusselt number. 

k. Calculate fin efficiency using previously- 
calculated, heat- transfer coefficient. 

l. Recalculate heat-transfer coefficient. 
Recalculate Nusselt number and compare to Nusselt 

number found in j and iterate if not within 1 % . Figure 5.1 
shows an example of plotted data and Appendix G shows a 
listing of PLOT. 

2. Program NSFDRP 

a. Compute average cooling water temperature. 

b. Compute cooling water velocity. 

c. Compute mass flow rate of cooling water. 

d. Compute heat transferred to cooling water. 

e. Compute log mean temperature difference. 

f. Compute overall heat- transfer coefficient based 
on outside area. 
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g. Compute wall resistance based on outside area. 

h. Compute Reynolds number of cooling water. 

i. Assume fin efficiency. 

j. Compute inside heat- transfer coefficient. 

k. Calculate cooling water temperature rise. 

l. Calculate fin efficiency using previously- 
calculated, heat- transfer coefficient. 

m. Recalculate heat-transf er coefficient and iterate 
if not within 1 % of the value found in h. 

n. Calculate condensing heat-transfer coefficient 
from Uq, R w and h^ . 



TABLE I 

HP 3497A Channel Assignments 



Channel 


Assignment 


20 


steam 


21 


steam 


22 


condensate 


23 


room 


24 


room 


25 


wall 


26 


wall 


27 


wall 


28 


wall 


29 


wall 


30 


wall 
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Figure 4«1 Photograph of Data Acquisition/Reduction System 








V. RESULTS AND DISCUSSIONS 



A. INSIDE HEAT-TRANSFER COEFFICIENT 

Experimental data were taken at atmospheric pressure by 
the procedures outlined in Chapter III. Twelve data sets 
were taken to determine the inside heat- transfer coefficient. 
The first four data sets were utilized to test the data 
acquisition/reduction program SIEDER and the middle five 
data sets experienced mixed filwise and dropwise condensa- 
tion conditions. Data sets 10-12 showed filmwise condensa- 
tion throughout each of the data collection sets. Figure 5*1 
is a plot of data taken during run number 10. Similar results 
were obtained for data runs 11-12. The analysis of all 
data showed a Sieder-Tate coefficient of 0.029-0.001 to be 
representative for the experimental apparatus (Table II). 

The Sieder-Tate correlation for fully developed turbulent 
flow has a leading coefficient of 0.027 {Ref. 1 6} . Rose 
{Ref. 14) shows a leading coefficient of 0.030 with an 
L/D ratio of 1 7 , which is in good agreement with the results 
of this thesis. The L/D ratio for the experimental apparatus 
used during this thesis was 18, well below that for fully 
developed flow {Ref. 1 6} ; therefore the inside heat- transfer 
coefficient was correspondingly higher. 
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TABLE II 



Sieder-Tate Coefficients 



File 


Ci 


1 0 


.0291 


1 1 


.0287 


12 


.0287 



B. OUTSIDE HEAT-TRANSFER COEFFICIENT 

Experimental data were taken at atmospheric pressure by 
the procedures outlined in Chapter III. Four data sets wer 
recorded to detemine the outside heat- transfer coefficient 
using the program NSFDRP (Appendix G) to reduce the raw 
data. Figure 5.2 shows a plot of the four data sets along 
with a curve showing the corresponding Nusselt prediction. 
All of the data lie above the theoretical curve. This 
discrepancy has two possible explanations. Either dropwise 
condensation and/or vapor shear associated with a vapor 
velocity of 3 m/s could increase the condensation heat 
transfer coefficient. The effect of vapor shear has been 
studied by many investigators including Rose as described 
in a very recent report {Ref. 18}. Rose reported that when 
vapor velocity was increased to 7 m/s for refrigerant-113, 
the condensing coefficient increased by 20%. The presence 
of dropwise condensation significantly raises the heat- 
transfer coefficient, sometimes as much as ten times 
{Ref. 17}. 
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A large amount of scatter occurred in the recorded data. 
This situation could have occurred due to noncondensable 
gases and the inability to regulate system pressure. During 
normal operation, the Mechanical Engineering Department 
air compressor ( 1 60 psi) was used to operate the air 
ejector. During the condensation heat transfer coefficient 
runs, this air compressor was inoperative; therefore, house 
air (50 psi) was used. This source of air had a very low 
capacity and fluctuated greatly during each data run, 
perhaps causing varying amounts of air to be in the vicinity 
of the test tube. 

Due to limited and lack of reliable data, very little 
can be concluded for the results thus far for the condensing 
coefficient. More data is necessary to reach any final 
conclusions . 
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Figure 5.1 Sieder-Tate Paramater Plot 
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Figure 5.2 Condensing Coefficient Plot 



VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1 . An experimental apparatus that allows the study of 
condensation heat transfer of steam on a single, horizontal 
tube has been designed, constructed and instrumented. 

2. The water-side heat-transf er coefficient was 
established using the Sieder-Tate correlation with a leading 
coefficient of 0.029 - 0.001. 

3. An analysis scheme to establish the condensing heat- 
transfer coefficient has been outlined, and data acquisition/ 
reduction programs have been written and tested. 

B. RECOMMENDATIONS 

1. To prevent leaking, install Swagelok fittings with 
Teflon ferrules into the dump condenser lines penetrating 
the stainless base of the dump condenser section. 

2. Rig suspension cables to remove the dump condenser 
section by attaching them in the ceiling using metal studs. 

3. Install a thermopile at the outlet of the test 

condenser tube as a second means of measuring the cooling 

water temperature rise across the test section. 

/ 

4. Use a commutating bridge to accurately determine 
boiler heater resistance over the full range of power 
settings . 
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5. Install a nitrogen cold trap on the mechanical vacuum 
pump. 

6. Install a permanent drain line on the cooling water 
reservoir to replace the temporary siphon tube. 

7. Replace the extension wire from the ice bath to the 
data acquisition system with thermocouple-quality copper 
wire . 

8. After recommendation No. 1 above is completed, check 
the system for leak tightness under vacuum conditions. 

9. Complete experimentation to determine the condensing 
heat-transfer coefficient at 1 atmosphere as well as under 
vacuum conditions. 
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APPENDIX A 



SAMPLE CALCULATIONS 

1 . INSIDE HEAT-TRANSFER COEFFICIENT 

A set of calculations was performed to verify data reduc- 
tion program SIEDER. One data point was selected arbitrarily 
from data set 10. 

a. Given the following inlet and outlet temperatures, 
compute average cooling water temperature. 

T 1 = 23 . 94°C 
T 2 = 26.32°C 

T, = (23.94 + 26.32)/2 = 25.13°C 
b 



b. Compute average wall temperature. 

T = (T ,+T 0 +T q +T , + T -+T ,)/6 = 62.64°C 
w wl w2 w3 w4 w5 wo" 

c. Compute water velocity.' 

Flow meter reading = 25% and from the calibration 

curve, V = 4.5 ft/s or V = 1.37 m/s. 
w w 

From Reference 1 6 , Table A. 6 at 300 K the following 
properties are found: 

?r = 5.83 

U f = 855E-6 N-s/m 2 
C^ = 4.179 kJ/Kg-K 
X f = 613E-3 W/m-K 
v f = 1 .003E-3 m 3 /kg 

d. Compute mass flow rate of cooling water. 

A = mD i 2 /4 

m = V t _ A/v = 0.173 kg/s 
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e. Compute heat transferred to cooling water. 
Q = m C p (T r T 2 ) 

Q = 1721 W 



f. Compute average inside wall temperature. 

T = Q In (D /D.) /2 it k. L Reference 17 

W O 1 X 

T = 2 . 1 6°C 
w 

T . T - 0. 5 AT 
Wl = w w 

T . = 61 . 56°c 

wi 

g. Compute log mean temperature difference. 

LMTD = (T 2 - Ti )/l n [(T w± - Ti >/(T wi -T 2 )j 

h. A fin efficiency is assumed to be 20% for first 
interaction. 

i. Compute heat- transfer coefficient, 

hn = Q/tt D. (1+^ +L 2 *n 2 ) LMTD 

h. = 7772 W/m 2 -K 

l 

j. Compute Nusselt number. 

Nu = h. D. /ic 

l i 

Nu = 1 61 first iteration 

k. Calculate fin efficiency. 

?, = 2tt D = 0.1 197 m 

1 o 

P 0 = 2tt D = 0.1 00 m 

2 r 

A, = tt D (D -D. ) = 0.00038 m 2 

1 OOl 

A 0 = it D (D -D. ) = 1 .6 E-4 m 2 

2 r r l 



= / h ^ P -] / k £>A ^ 



m 1 = 79.74 

"2 = /h i P 2 /k f A 2 
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m 2 = 112.32 

H 1 = tanh (m 1 )/m L 

n 1 = -2079 
n 2 = .2913 

l. Recalculate heat-transf er coefficient. 

= 7589.9 W/m 2 -K 

m. Recalculate Nusselt number and compare to previous 
calculated value. 

Nu = 157.24 

This value for the Nusselt number agrees within 0.3% of 
the computer generated value (Nu^ = 156.7). 

2. CONDENSATION HEAT-TRANSFER COEFFICIENT 

A set of calculations was performed to verify the data 
reduction program NSFDRP. One data point was selected from 
data run 1 . 

a. Given the following inlet and outlet cooling water 
temperatures, compute average water temperature. 

T 1 = 30 . 06°C 

T 2 = 32.53°C 

T b = (T.j +T 2 )/2 = 31 • 3°C 

T = 96.7°C 
s 

b. Compute water velocity. 

Flow meter reading = 30% and from the calibration curve 

V = 1.64 m / s 

w 

c. Compute mass flow rate of cooling water. 

From Reference 16, Table A. 6 at 304& the following 
properties are found: 
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Pr = 5.29 

U = 0.781 E-3 N-s/ra 2 
C p = 4.178 kJ/kg-K 

k = 0.62 W/m-K vw = 1.005 E-3 m 3 /kg 
A = ttD. 2 /4 

1 ' 

A = 1.267 E-4 m 2 

m = V A/ v 
w 

m = .2067 kg/s 

d. Compute heat transferred to cooling water. 

Q = 2133 W 

e. Compute LMTD. 

LMTD = (12-T^/ln [(T s -T 1 )/(T s -T 2 )J 
LMTD = 65.40°C 

f. Compute overall heat- transfer coefficient based on 

outside area. 

A = 7T D L 
o o 

U q = Q/A c LMTD 

U = 4086.71 W/ra 2 -K 
o 

g. Compute wall resistance based on outside area. 

R = D In (D /D . ) /2 k , 
m o o i cu 

R = 1 .003 3-5 m 2 -K/W 
m 

h. Compute Reynolds number of cooling water. 

Re = V D ./v, 
w w 1 f 

Re = 26535 
w 

i. Assume fin efficiencies of = 0.20 and n 2 = 0.24 
and assuming (u/u w ) ^ ^=1 « 0 
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j • 



h. 



1 

h. 



1 



Compute inside heat- transfer 
= (k w C./D.) Re u - 8 Pr' 333 (u/u 
= 8533.72 W/ra 2 -K 



coefficient. 




U 



k. Calculate fin efficiency using previously-calculated 
heat- transfer coefficient from j. 
m^ = 83 . 56 m""' 

m 2 = "^7.70 m 

n 1 = 0.20 

n 2 = 0.24 



i. 



T 

w 



Calculate cooling water rise. 

= T. + AT = 72° C 
b w 



From Reference 1 6 at 345 K, v = 389E-6 N-s/m' 

w 



(u/ u M ) 



0.14 



= 1.1 

m. Recalculate heat transfer coefficient and iterate if 
not within 1 % of the value found in j. 
h. = 9387 W/m^-K 



n. Calculate condensation heat- transfer coefficient. 

1 



h 



1 



D L 
o 



U 



D i (L+n 1 L 1 +n 2 L 2' )h i 



m 



h = 1 0405 W/m 2 -K 
o 

This value for the condensation heat-transf er coefficient 
agrees within 1 % of the computer generated value (h Qc = 1 0550). 
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APPENDIX B 



SYSTEM STARTUP PROCEDURES 

1 . Familiarize all valve numbers by referring to Figure 2.1 . 

2. Fill boiler with distilled water to a level 4 inches 
below top of sight glass. 

a. Start vacuum pump or air ejector. 

b. Open suction line supply valve (3). 

c. Allow vacuum pump/air ejector to run until 15 inches 
Hg of vacuum is attained. 

d. Close suction line supply valve (3). 

e. Shut down vacuum system. 

f. Open discharge on distilled water tank. 

g. Open distilled water supply line valve (6). 

h. Close distilled water supply line valve (6) when 
boiler is filled to 4 inches below the top of the 
sight glass. 

3. For vacuum operations, start vacuum pump/air ejector 
then open suction line supply valve (3) until desired 
pressure is reached. Suction line supply valve may be 
throttled to maintain pressure. 

4. Energize viewing window heater and adjust rheostat to 
110. Open air supply valve (9) about g turn. The air 
supply may require further adjustments to keep viewing 
window clear. 

5. Open dump condenser supply valve (2) fully. 
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6. Fill cooling water reservoir. 

a. Open reservoir stop valve (10). 

b. Open test section throttle valve (11). 

c. Open reservoir fill valve and allow tank to fill to 
about 101 .6 mm (4 in) from the top. 

d. Close reservoir stop valve. 

e. Energize cooling water circulation pump. The switch 
is located above the pump. 

7. Shut the following power supply breakers: 

a. Panel p-5 switch 3. 

b. Switch 3 on side of control panel. 

c. Switch 1 on front panel. 

8. Increase boiler- heater voltage to a maximum of 200 
volts during warmup; about 15 minutes. When operating under 
vacuum conditions use a lower voltage to limit vibration 

of glass piping. 

9. As steam is generated, air is expelled through a relief 
valve. System pressure should never be allowed to exceed 
10 psig. 

10. Once steam reaches condenser test section (about 30 
minutes following startup), full power operation may begin. 

11. The apparatus should be allowed to steam for 30 minutes 
prior to taking any data. 
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APPENDIX C 



SYSTEM SHUTDOWN PROCEDURES 

1. Decrease boiler heater voltage to zero. 

2. Open all three circuit breakers. 

3. Secure view window heater rheostat switch; allow air 
to flow for about 15 minutes longer. 

4. Secure cooling water circulation pump. 

5. Close cooling water throttle valve. 

6. Close dump condenser supply valve. 

7. If in use, secure vacuum suction supply valve followed 
by vacuum pump/air ejector. 

8. For cleaning purposes only, the boiler may be drained 
by unscrewing the Swaglok nut housing the condensate return 
thermocouple. The Teflon ferrules require only hand 
tightening. 
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APPENDIX D 



TEST SECTION INSTALLATION PROCEDURES 

1. Install outlet nylon holder, mixing chamber and return 
line onto provide 6.35 mm (0.25 in) diameter bolts and 
tighten . 

2. Place clean test tube into inlet nylon holder and care- 
fully slide holder and test tube into condenser test section 
snugly. 

3. Align test tube to have one thermocouple channel vertical. 

4. Carefully slide thermocouples into channels. 

5. Ensure wires are straight and lined up with markings on 
inlet nylon holder and tape into place. 

6. Position nylon thermocouple wire retainer over wires. 

7. Place inlet nylon holder over end of test tube and 
tighten all bolts. 

8. Energize circulation pump and check for leakage around 
connections both inside and outside test section. 
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APPENDIX E 



THERMOCOUPLE CALIBRATION 

1 . EQUIPMENT USED 

a . Thermocouple wire 

Copper-constantan , 0.245 mm (0.01 in) Teflon coated wire 
was used for all thermocouples. (Omega Catalogue) 

b. Calibration bath 

A Rosemount Engineering Co. Model 91 3A calibration bath 
shown schematically in Figure A.1 was used. 

1) Heating: Electrical 

2) Cooling: Liquid Nitrogen 

Note: 

Once a desired temperature is reached, the temperature 
is held constant by rapid cycling between heating and cooling. 
The manufacturer rates the temperature bath fluctuation 
to be less than 0.01°C. Observations during the calibra- 
tion procedure showed bath temperature fluctuation to be 
less than 0.02 5°C (l.OqV). 

c . Thermocouple readout 

An HP 3054A Automatic Data Acquisition/Control System 
with resolution of 0.1 uV was used to obtain data. 

d. Reference Point 

A reference junction was placed in a 0.15 m Dia X 0.28 m 
high Thermos flask filled with finely-crushed ice. 
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e. Temperature measurement 

Mercury-in-glass thermometers with a resolution of 0.1°C 
were used. A total of four thermometers were used with the 
following ranges: 0-25, 25-54, 50-76 and 76-1 04°C. 

2. PREPARATION 

a . Procedure 

Two thermocouples were prepared, one at the beginning 
and the other at the end of a reel of thermocouple wire, 
using the following procedure: 

1 ) The Teflon insulation was removed for a length of 
about 4 ram from one end of a 2 m long piece of wire 
and a thermocouple bead was made using a Dynatech 
Corporation thermocouple welder. 

2) The other end of the thermocouple wire was soft 
soldered to a pair of heavy duty copper wires. 

3) Heat shrink insulation was placed on the soldered 
end. 

b. Reference point 

The thermocouple reference junctions were placed in a 
5 cm test tube and then into a thermos flask containing 
finely crushed ice as shown in Figure A.1. 

c . Thermometers 

Thermometers were completely immersed in the calibration 
bath as recommended by the manufacturer. 

d. Analysis 

A short computer program (TC_CAL) was written to accept 
the thermocouple readings through the Hp data acquisition 
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system, and the bath temperature through the keyboard. This 
program prints all data as well as stores the data on a 
computer disk. A listing of the computer program is provided 
in Appendix G. 

3. CALIBRATION PROCEDURE 

a. The bath temperature was set at about 15°C. 

b. A constant temperature is reached when rapid cycling 
of heating and cooling occurs, but wait three minutes to 
enter data. 

c. Enter the bath temperatures, as measured by the 
mercury-in-glass thermometer, into the keyboard. Immediately 
following the data entry, the computer will read the 
thermocouple EMF automatically and print the bath temperature 
and thermocouple EMF (uV). 

d. When the computer prompts for more data, increase 
bath temperature by about 10°C and repeat the above procedure 
up to about 95°C. 

e. Repeat the above procedure for decreasing temperatures. 

4. CALIBRATION CURVE 

a. A least-squares technique was used to generate a 

third order polynomial in the form of: 

2 3 

= a^ji+a^i +a^n 

T is the temperature in degrees Celsius and E is the 
thermocouple reading in microvolts. 

b. Values generated for the coefficients are: 
a., = 0.02635 
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a 2 = -9.735 E-7 
a 3 = 6.577 E-11 

c. Figure A. 2 is a sample printout of calibration 
data . 

d. Figure A. 3 is a graphical representation of calibra- 
tion data. 
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Figure E.1 Schematic of Calibration Bath 



Month, date and t ime : 03 : 1 6 : 20 : 20 : 30 



Bath Temp 
(Deg C) 

13.06 

24.33 

35 . 33 

46.72 
59.50 

72.06 

34.73 

94.00 
33 . 22 
82.89 
66.39 

54.33 

42.06 
30.67 
20.72 

13.06 



16 runs 



Enf (0) 
(Micro V ) 

503.0 

955.0 

1409.5 

1839.8 

2442.9 

2996 . 4 

3575.3 

3998.8 

3963.9 

3483.3 

2767.8 

2216.9 
1691 .7 

1213.4 

806.5 

502.5 



c n r ( ! ) 
(Micro V > 
503.2 

956.0 

1409.0 

1391 .1 
2*42.4 

2995.7 

3575. 1 

3993.1 

3963.0 

3433.2 

2766.9 

2215.8 

1690.9 

1212.4 

806.0 

502.0 



were stored in file C£L_DATA 



Figure E.2 Calibration Printout 
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Figure E.3 Thermocouple Calibration Curve. 



APPENDIX F 



TEST TUBE CLEANING PROCEDURES 
Caution 

The following mixture is highly toxic and very irritant 
when contacted by the skin. Be very careful not to splash 
the solution in any way. Wear protective gloves and eye 
protection when the actual cleaning is being conducted. 

1 . Ensure the surface of the tube is smooth by placing 
in lathe and polishing with copper cleaner (Brasso) and a 
clean rag. 

2. While still in the lathe, wipe completely with 
acetone followed by ethyl alcohol. 

3. Prepare a solution of equal parts of 50% sodium 
hydroxide solution (commercially prepared, available through 
MCB Manufacturing Chemists Inc., Cincinnati, OH) and 

ethyl alcohol. About 4 ounces each mixed in a shallow 
stainless-steel pan is sufficient. The sodium hydroxide 
will precipitate out of solution and form a white paste. 

4. Heat the mixture to about 80°C, but try to avoid 
spattering. The white pasty texture os the solution should 
now dissolve into a clear liquid. 

5. Completely immerse the tube into the mixture and 
scrub the entire surface with a bristle brush (old tooth- 
brushes work fine). 

6. Rinse with tap water followed by distilled water. 



61 



7. Allow the tube to remain in distilled water until 



ready for installation. 

8. The nylon holders should also be cleaned with the 
mixture . 

9. The test tube should be installed as soon as possible 
after cleaning. 
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APPENDIX G 



COMPUTER PROGRAM LISTINGS 



100! FILE NAME : TC_CAL 
110! REVISED: Jane 1. 1983 

120! 

130 DIM Emf(l) 

1 40 PRINTER IS 701 
! 50 BEEP 

160 INPUT "ENTER MONTH. DATE AND TIME (MM :DD:HH :MM:SS>" , Date$ 

170 OUTPUT 709 : "TD" : Dates 
ISO OUTPUT 709 .-"TD" 

190 ENTER 709: Dates 

200 PRINT USING " 1 2X . ""Month . date and time : "" , 1 4A" :DateS 
210 BEEP 

220 INPUT "GIVE A NAME FOR DATA FILE" . D_f i leS 

230 CREATE BOAT D_file$.5 

240 ASSIGN s>F i 1 e TO D_fileS 

250 J=0 

260 PRINT " " 

270 PRINT USING "!2X.""Bath Temp Emf<0) Emfd) 

280 PRINT USING "12X."" (Deg C) (Micro V) (Micro V) 

290 Repeat: ! 

300 BEEP 

310 INPUT "ENTER BATH TEMPERATURE (DEG F)",T_bath 
320 T_bath= < T_bath-32 . ) /I . 8 
330 OUTPUT 709 : "AR AF20 AL21" 

340 OUTPUT 722: "FI R1 T1 Z1 FL1 " 

350 FOR 1=0 TO I 
36C OUTPUT 709: "AS SA" 

370 ENTER 722:Emf(I> 

380 Emf (I)=ABS(Emf ( I ) ) *1 0~S 
390 NEXT I 

400 PRINT USING " 1 3X . 3D . DD . 7X . 4D . D . 5X , 4D . D" : T_bath . Emf ( *) 

410 OUTPUT i«Fi le ; T_bath ,Emf ( *) 

420 BEEP 

430 INPUT "ARE YOU TAKING MORE DATA ( 1 =YES . 0=NO>?" .Go_on 
440 J=J+1 

450 IF Go_on=1 THEN Repeat 

4bU BtLP 

470 PRINT " " 

480 PRINT USING "12X.DD."" runs were stored in file "" . 1 0A" : J , D_f i le'S 
490 END 
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<000! FILE NAME: NSFDRP 
1010 ! REVISED: 23 MAY 1983 

1020 DIM Emf (10) 

1030 Di- .01 27 
1040 Do-. 01 905 
1050 Dr-. 015875 
1060 Dssp= . 1 524 

1070 Ax = PI*Dssp '2/4-PI*Do*L ! TO BE MODIFIED 

1080 L=. 13335 

1090 LI -.060325 

1100 L2-. 034925 

1110 L 1 =L 1 +< Do-D i >/4 

1120 L2-L2+ ( Dr-D i ) /4 

H30 Ci=.029 

1140 Kcu-385 . 

1150 Rm=Do«L0G(Do/Di>/(2*Kcu> 

1160 PRINTER IS 701 
1170 CLEAR 709 
H80 BEEP 

1190 INPUT "ENTER MONTH. DATE AND TIME (MM : DD : HH :MM : SS) ". DateS 
1200 OUTPUT 709 : "TD" : DateS 
1210 OUTPUT 709 : "TD" 

1220 ENTER 709;Date$ 

1230 PRINT " Month, date and time :";Date$ 

1240 BEEP 

1250 INPUT "ENTER INPUT MODE ( 1 -3054A .2-FILE)" . Im 
1260 IF Im-1 THEN 
1270 BEEP 

1280 INPUT "GIVE A NAME FOR THE RAW DATA FILE" . D_f i leS 
1290 CREATE BDAT D_fileS,10 
1300 ELSE 
1310 BEEP 

1320 INPUT "GIVE THE NAME OF THE EXISTING DATA FILE" . D_f i leS 
1330 PRINT " " 

1340 PRINT USING "!2X.""7he following analysis was performed for data in file 
" . 10A":D_f ileS 
1350 BEEP 

1360 INPUT "ENTER THE NUMBER OF RUNS STORED". Nrun 
1370 END IF 
1380 BEEP 

1390 INPUT "GIVE A NAME FOR PLOT DATA FILE" ,P_f i leS 

1400 CREATE BDAT P_fileS.5 

1410 ASSIGN SFile TO D_fileS 

1420 ASSIGN SFilep TO P_fileS 

1430 J-0 

1 440 Repeat : ! 

1450 J-J+1 

1460 IF Im- 1 THEN 

1470 BEEP 

1480 INPUT "ENTER FLOWMETER READING”. Fm 
1490 OUTPUT 709 : "AR AF20 AL30" 

1500 OUTPUT 722 : "F ! R1 T1 Z1 FL1" 

1510 OUTPUT 709: "AS SA” 

1520 ENTER 722:B P v 
1530 FOR 1=0 TO 10 
1540 OUTPUT 709: "AS SA" 

1550 IF I>4 THEN 

1560 Se-0 

1570 FOR K-0 TO 10 
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1580 


ENTER 722 ;E 




1590 


Se=Se+E 




1600 


NEXT K 




1610 


Emf(I)=ABS(Se/I0>*1 .E+6 




1620 


ELSE 




1630 


ENTER 722 : E 




1640 


Emf ( I>=ABS(E>*1 .E+6 




1650 


END IF 




1660 


NEXT I 




1670 


ENTER 722 : Emf ( I ) 




1680 


OUTPUT 7 1 3 : "T1 R2E" 




1690 


WAIT 2 




1700 


ENTER 7 1 3 ; T 1 




1710 


OUTPUT 713: "T2R2E" 




1720 


WAIT 2 




1730 


ENTER 7 1 3 ; T2 




1740 


IF J=1 THEN 




1750 


BEEP 




1760 


INPUT -ENTER NANOMETER READING (mm Hg)",Phg 




1770 


BEEP 




1730 


INPUT "ENTER HEIGHT OF MANOMETER WATER COLUMN" . Pmater 




1790 


END IF 




1800 


ELSE 




1810 


ENTER 9File:Emf(*),Fm,T1,T2,Phg. Puater , Bpv 




1820 


END IF 




1830 


Is teamed 




1840 


FOR 1=0 TO 1 




1850 


Tsteam=Tsteam+ . 5*FNTvsv(Emf ( I) ) 




1860 


NEXT I 




1870 


Tum^O * 




1880 


FOR 1-0 TO 5 




1890 


T iji < I)=FNTvsv<Emf(I+5)> 




1900 


Tum=Tidm + Tu( I ) 




1910 


NEXT I 




1920 


Tum=Tum/6 




1930 


Tcon=FMTvsv(Emf (10)) 




1940 


Psat=FNPvst ( Ts team) 




1950 


P tes t * < Phg+Puater/ 1 3 . 6 ) * 1 33 . 322 




1960 


IF J=1 THEN 




1970 


Vst=FNVvst < Ts team) 




1980 


Ppng=( Ptest-Psat ) /Ptest 




1990 


Ppst= 1 -Ppng 




2000 


f1tng = Ppng/<287*<Tsteam^273« J 0 ) / C Vs t *P tes t > +Ppng ) 




2010 


Vfng=Mfng/(1 , 808- . 608*Mf ng) 




2020 


Mf ng = Mf ng* 1 00 




2030 


Vf ng = Vf ng*1 00 




2040 


BEEP 




2050 


PRINT M M 




2060 


PRINT USING *' 1 2X , "’’Measured pressure 


= "" . D . 3D 


r* «« 1 1 

c. . 


(Pa) ’:Ptest 




2070 


PRINT USING " 1 2X , ""Pressure corresponding to measured steam 


temp = "" , D . 3D 


E. 


(Pa) :Psat 




2080 


PRINT USING " 1 2X , ""NoncondensaJble gas concentration 


= "" . D . D , "" ( 


Vol 7) : Vf ng 




2090 


PRINT USING "12X. ""Noncondensable gas concentration 


= ”" . DD . D , ""( 


Mass 


7) :Mfng 




2100 


IF Mf ng> . 5 THEN 




21 1 0 


BEEP 




2120 


PRINT " " 




2130 


IF Im=1 THEN 
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2 MO BEEP 
2150 PRINT " " 

2160 PRINT USING " 1 OX , ‘"'Energi ze the vacuum system 
2170 BEEP 

2130 INPUT "OK to accept the present run ( 1 = YES , G=NO> ?" . Ok 
2130 IF 0k=0 THEN 
2200 BEEP 

2210 DISP "NOTE: THE PRESENT DATA RUN WILL BE DISCARDED!! " 

2220 WAIT 5 

2230 GOTO M60 

2240 END IF 

2250 END IF 

2260 END IF 

2270 END IF 

2280 IF Im= ! THEN OUTPUT 3Fi le :Emf<*> ,Fm. T! , T2 .Phq .Pwater .B ps 
2230! ANALYSIS BEGINS 
2300 Tavg=< Tt +T2)*.5 
2310 Cpu=FNCpu(Tavg) 

2320 Rhoij=FNRhou( Tavg) 

2330 Mf=Fm*. 02451*. 3048~3/100 

2340 Md=Mf *Rhou 

2350 Vu-Nf/(PI*Di“2/4) 

2360 Q = Md*Cp(j*< T2-T1 ) 

2370 Qp=0/(PI*Do*L) 

2380 Tijo = Twm+Qp*Rm* .5 
2330 Dtf =Tsteam-Tuo 
2400 Ku=FNKu( Tavg) 

2410 Muu=FNMuu( Tavg> 

2420 Reu)=Rhoij*Vu*Di/Muu 
2430 Pru=FNPru( Tavg) 

2440 Fe 1 =0 . 

2450 Fe2=0 . 

2460 Of = t . 

2470 Hi=Ku*Ci/Di*Rew ' .3*Pru' .3333*Cf 
2430 Dt=G/<PI*Di-<L+L1*Fe1 +L2*Fe2)*Hi> 

2430 Cf c= (Nuw/FNMuu( Tavg + D t > ) ' . 14 

2500 IF ABS( (Cfc-Cf )/Cfc>> .01 THEN 

2510 Cf-<Cf+Cfc>*.5 

2520 GOTO 2470 

2530 END IF 

2540 P1=PI*(Di+Do) 

2550 A1=<Do-Di)*PI*(Di+Do)*.5 

2560 til = (Hi*Pl /(Kcu*A1 ) )" .5 
2570 P2=PI*( Di+Dr) 

2530 A2=(Dr-Di)*PI*(Di+Dr>*.5 

2530 M2= < Hi*P2/ (Kcu*A2> > ' .5 

2600 Lmtd=(T2-Tl )/L0G< ( Ts team-TI > / < Ts team- T2) ) 

2610 Uo=0/(Lmtd*PI*Do*L) 

2620 Fe1=FNTanh(M1*L1)/(M1*L1> 

2630 Fe2=FNTanh(M2*L2)/<H2*L2) 

2640 Ho= 1/(1 /Uo-Do*L/(Di*<L+Ll *Fe1 +L2*Fe2) *Hi ) -Rm) 

2650 Dtc=Q/<PI*Di*(L+L!*Fe1 +L2*Fe2)*Hi > 

2660 IF ABS( (Dtc-Dt)/Dtc)> .01 THEN 2470 
2670 Troom-FNTvsvUEmf (3>+Emf (4> )*.5) 

2680 PRINT "U0=" :Uo 

2630 Hos=0/( p I*Do*L*< Tsteam- Two) ) 

2700 PRINT "H0_STAR=" ;Hos 
2710 Hf g=FNHf g < Ts team ) 

2720 Tf i lm = ( Ts team* I uo) * . 5 
2730 Kf-FNKu(Tfilm) 
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2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
Hqp 
2860 
2870 
2880 
2880 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
D_f ileS 
2970 END IF 
2980 BEEP 
2S90 PRINT " " 

3000 PRINT USING M 1 2X , ""NOTE : , DD , ”"X- Y pairs were stored in plot data file M 

".10A":J.P_f ileS 

3010 ASSIGN ?File TO * 

3020 END 

3030 DEF FNPvstf Tsteajn) 

3040 DIM K ( 8 ) 

3050 DATA -7 . 69 1 234564 . -26 . 08023696 1 58 . 1 706546 . 64 . 23285504 1 1 3 . 9646225 
3060 DATA 4 . 1 67 1 1 732 , 20 . 9750676 , 1 . E9 ,6 
3070 READ K<*) 

3030 T=<Tsteam+273. 15)/647 .3 

3090 Smn = 0 

3100 FOR M=0 TO 4 

3110 Sum=Sum+K(N)*( 1-D '(N+1) 

3120 NEXT N 

3130 dr = Suxi/(T*(1-rK(5)*(1-T) : K(6 » *< 1 -T> '2) >-< 1 -T) / (K( 7 )*< i -T ) '2+K(8> > 

3140 Pr=EXP(Br ) 

3150 P=22 1 20000*Pr 

3160 RETURN P 
3170 FNEND 
3180 DEF FNHfg(T) 

3190 Hfq=2497.7389-T*(2.2074+T*U . 7079E-3-2 . 8593E-6 ) ) 

3200 RETURN Hfg*1000 

3210 FNEND 

3220 DEF FNMuu(T) 

3230 Mu-! .5760947 3E-3-T*<3. 51 1 98576E-5-T*<3 . 58358 1 6E-7- 1 . 365536 1 1 5E-9*T> > 

3240 RETURN Mu 

3250 FNEND 

3260 DEF FNVvst(T) 

3270 0=58 . 4525538 -T*( 1 .5 1 508776-W . 0 1 372746585-T*4 .253667 1 1 E-5 ) ) 

3280 RETURN 0 
3290 FNEND 
3300 DEF FNCpu(T) 



Rhof=FNRhow(Tf ilm) 

Muf =FNMuu< T film) 

Hnu= . 725* < Kf "3*Rhof '2*9 . 8 1 *Hf g/ < Muf *Do*D tf ) > ~ . 25 
Hpg= .65l*Kf*(Rhof~2*9.81*Hfg/( Muf *Do*Qp ) ) ^ .3333 



PRINT 

PRINT 

PRINT 

PRINT 



USING 

USING 

USING 



PRINT USING 
PRINT USING 



PRINT 

PRINT 



USING 

USING 



"12X. ""Water velocity 
"1 2X . ""Reyno Ids number 
" I2X , ""Inside heat-tran coef 
" 1 2X . ""Heat flux 
" 1 2X . ""Outs ide heat-tran coef 
1 2X . ""Nusse 1 t coefficient 



" . Z . DD , <m/S> ; Vu 

" . 5D . D" : Reu 

" , 5D . D , (H/m'2-K) ; Hi 

" . Z . 5DE . "" <W/m"2> :Qp 

" . 5D . D . "" (W/m'2-K) ;Ho 

. 5D . D . "" <M/n'2-K) ;Hnu 



" 1 2X . 1 2X , ""Nus coef (from heat flux) = "" . 5D>D , "" < W/m '2-K )' 



OUTPUT ?F i 1 ep : Op , Hos , Hnu 

IF Im-1 THEN 

BEEP 

INPUT "ARE YOU TAKING MORE DATA ( 1 =YES . 0=110) ?" ,Go_on 

IF Go_on=1 THEN Repeat 

ELSE 

IF J<Nrun THEN Repeat 
END IF 

IF Im-1 THEN 
BEEP 

PRINT USING " 1 2X , ""NOTE : ,,M ,DD,"" data runs were stored in file "".I0A";J, 
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3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3330 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 



Cpu=4 .2 1 1 20858-T*<2 . 26826E-3-T*<4 . 4236 1 E -5+2 . 7 1 428E-7*T> > 

RETURN Cpu*1 000 
FNEND 

DEF FNRhou(T) 

Ro = 999 . 52946+ T*< . 0 1 269-T»<5 . 4825 1 3E - 3 - T-* t .234147E-5)) 

RETURN Ro 
FNEND 

DEF FNPru(T) 

Prw= 1 0~ < 1 . 09976605-T*< 1 .3749326E-2-7*<3.968875E-5-3.45026E-7*T>>> 

RETURN Prut 

FNEND 

DEF FNKu(T) 

Ku= . 5625894 +T*< 2. 2964 546E-3-T*< 1 .509766E-5-4 .0581652E-8*T) ) 

RETURN Ku 
FNEND 

DEF FNTanh(X) 

P=EXP(X) 

Q=EXP(-X> 

Tanh=(P+Q)/(P-Q> 

RETURN Tanh 
FNEND 

DEF FNTvsv(V) 

T=V»( . 026 1 73344 1 6- V* (9 . 2447859E-7-V*6. 0746642E- 1 1 ) > 

RETURN T 
FNEND 

DEF FNHf(T) 

Hf=T*(4.203849-T*(5.88132E-4-T*4.55160317E-6)> 

RETURN Hf * 1 000 
FNEND 
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1 000 ! 
1010 ! 
1020 
1030 
1040! 
1050! 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1 150 
1160 
1170 



FILE NAME: SIEDER 
REVISED: May 23. 1983 

DIM Emf < 1 0 > , Tu(5> 

Kcu=385 

Assign Geometric Variables 



Di-,0127 JInside Diameter 
Do=. 01905 !0utside Diameter 



Dr=. 015875 [Outlet End Diameter 
L= . 13335 !Condensing Length 
LI*. 060325 ! Inlet End "fin" Length 

L2=. 034925 !Outlet End "fin" Length 
PRINTER IS 701 
BEEP 

CLEAR 709 

INPUT "ENTER MONTH. DATE AND TIME (MM : DD : HH : MM :SS>" , 
OUTPUT 709;"TD":BS 
Series: ! 



B$ 



1180 OUTPUT 709;"TD" 

M90 ENTER 709 :A$ 

1200 PRINT USING " 1 2X . "".Month . date and time: "".14A";AS 
1210 BEEP 

1220 INPUT "ENTER INPUT MODE ( 1 =3054A ,2-FILE) ",1m 
1230 IF Im-1 THEN 
1240 BEEP 

1250 INPUT "PROVIDE NAME FOR THE DATA FILE" . D_f i leS 
1260 CREATE BDAT D_file$,10 
1270 ELSE 
1280 BEEP 

1290 INPUT "PROVIDE NAME OF THE DATA FILE" . D_f i le$ 

1300 BEEP 

1310 INPUT "ENTER THE NUMBER OF RUNS STORED". Nrun 

1320 PRINT USING "12X,”"The following analysis was performed 

". l0A":D_f lies 

1 33C END IF 

1340 BEEP 

1350 INPUT "PROVIDE NAME FOR PLOT DATA FILE". Plots 
1360 CREATE BDAT Plots, 5 
1370 ASSIGN SFile TO D_file$ 

1380 ASSIGN @Filep TO Plots 
1390 J=0 



i4fin Sxs = 0 
1410 Sxy =0 
1420 IF Im-1 THEN 

1430! READ DATA THROUGH DATA ACQUISITION SYSTEM 
1440! IF THE INPUT MODE (Im) = 1 
1450 BEEP 

1460 INPUT "ENTER FLOWMETER READING". Fm 
1470 OUTPUT 709 : "AR AF20 AL30" 

1480 OUTPUT 722 ; "FI R1 T1 Z1 FL 1 " 

1490 FOR 1=0 TO 10 
1500 OUTPUT 709: "AS SA" 

1510 IF I>4 THEN 



1520 Se=0 



for data in file 
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1530 FOR K=0 TO 9 

1540 ENTER 722 :E 

1550 Se=$e+E 

1560 NEXT K 

1570 Emf ( I )=ABS(Se/1 0) 

1580 ELSE 

1590 ENTER 722 :E 

1600 Emf < I ) =ABS(E> 

1610 END IF 
1620 NEXT I 

1630 OUTPUT 7 1 3 ; "T1 R2E" 

1640 WAIT 2 

1650 ENTER 7 1 3 ; T 1 

1660 OUTPUT 713; "T2R2E” 

1670 WAIT 2 
1630 ENTER 7 1 3 : T2 
1690 ELSE 

1700! READ DATA FROM USER-SPECIFIED FILE IF 
1710! INPUT MODE dm) =2 
1720 ENTER 9F i le ; Emf <* ) , T ! , T2 ,Fm 
1730 END IF 

1740! Compute Average Bulk Water Temperature 
1750 Tavg= ( T 1 +T2 ) * . 5 

1760! Compute Average Wall Temperature 

1770 Tuall-0 

1780 FOR 1=5 TO 10 

1 790 Ttiid-5) = FNTvsv(Emf (I)) 

1800 Twal l = Tuall+Tud-5> 

1810 NEXT I 

1820 Twal l=Twall/6 

1330! Compute Thermophysical Properties 
1840 Cpw=FNCpw( Tavg) 

1850 Rnow=r NRhow( Tavg) 

1860 Kw=FNKw( Tavg) 

1870 Mf = Fm*- . 0245 1 * . 3048 '3/ 1 00 !Volume Flowrate of Water (nT3/s) 

1880 Md=Mf*Rhou !Mass Flowrate of Water (kg/s) 

! '890 Vw=Mf /<PI*Di'2/4) !Water Velocity (n/s) 

1900 0=Md*Cpw*(T2-T1 ) !Heat Transferred to Cooling Water (W) 

1910 Dtw'Q^LGGtDo/DiJ/d’PI^Kcu^L) !Temp Drop Across Tube Wall (Deg C) 
1920 Twal l=Tual 1-Dtw* .5 

1930 Lmtd=(T2-T! ) /L0G< ( Twal 1-T I >/(Twall-T2) > 

1940 P1=PI»2*Do !Heat-Transf er Perimeter At Inlet 
1950 P2=PI*2*Dr ! Heat-Transf er Perimeter At Exit 
iSmO A1 = ( Do-D i >*PI*Do ! Heat Trar.sf c r Area At Inlet 
1970 A2= ( Dr-D i ) *P I*Dr ! Heat-Transfer Area At Exit 

1980 Hi=Q/(PI*Di*L*Lmtd) 

1990 Ml = (Hi*P 1 / (Kcu*A 1 ) ) " .5 
2000 M2=<Hi*P2/(Kcu*A2>)d5 

20'0 Fel =FNTanh(M1 *L 1 ) /( Ml *L 1 > !Fin Efficiency At Inlet 
2020 Fe2=FNTanh( M2*L2) / < M2*L2) !Fin Efficiency At Exit 
2030 Hic-Q/(PI*Di*(L+LI*Fe1+L2*Fe2)*Lmtd> 

2040 IF ABS( <Hi-Hic)/Hic)> .01 THEN 

2050 H i = (H ic+H i ) * . 5 

2050 GOTO 1990 

2070 END IF 

2030 PRINT " M 

2090 PRINT USING “ 12X . ""Pos i t ion number : V 2 3 4 

6 

2100 PRINT USING " 1 2X . ""Wal 1 temperature (Deg C) : ,6< DD . DD . 1X)";Tw(*) 

2110 PRINT USING "1 2X , "“Average wail temperature = **" , DD . DD . (Deg O"" 



5 

; Twal 1 
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2120 PRINT USING "? 2X , MM Inlet water temperature = "" . DD . DD . "" (Deg C) : T t 

2130 PRINT USING "12X , ’’"Out let water temperature = "".DD.DD,"" < Deg C)""":T2 

2140 PRINT USING " 1 2X . ""Log-mean- temp difference = "".DD.DD,"" (Deg C> :Lmtd 

2150! CALCULATE THE NUSSELT NUMBER 

2160 Nu=Hic*Di/Kw 

2170 Re=Rhow*Vw*D i /FNM udavg ) 

2180 Cf-<FNMu<Tavg>/FNMu<Tuall>>~. M 
2190 Prw=FNPrw( Tavg) 

2200 X=Re ' . 8*Prw ' . 3333*Cf tCompute Sieder-Tate Parameter 
2210! COMPUTE COEFFICIENTS FOR THE LEAST-SQUARES-FIT 
2220! STRAIGHT LINE 

2230 PRINT USING "12X, ""Water velocity = "".Z.DD,"" <m/S> ;Vw 

2240 PRINT USING " 1 2X , ""Sieder-Tate parameter = "",5D.D";X 

2250 PRINT USING " 1 2X . ""Nussel t number = "" . 4D . DD" ; Nu 

2260 OUTPUT ®Fiieo:X.Nu 

2270 Sx-Sx+X 

2230 Sy=Sy +Nu 

2290 $xs=$xs+X*X 

2300 Sxy =Sxy+X*Nu 

2310! STORE RAW DATA IN A USER-SPECIFIED FILE IF 
2320! INPUT MODE (Im) = 1 

2330 IF Im= 1 THEN OUTPUT «Fi lejEmt <*> , T1 ,T2,Fm 

2340 BEEP 

2350 J = J+1 

2360 IF In-1 THEN 

2370 INPUT "ARE YOU TAKING MORE DATA < 1 =YES , 0=NO>?" ,Go_on 
2380 Nrun=J 



2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2430 
2490 
2500 
2510 
2520 
2530 
n . D_f 

/. O 
L. U 

2550 

2560 

2570 

2530 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2630 

2690 

2700 



IF Go_on= 1 THEN 1420 
ELSE 

IF J < N r un THEN 1420 
END IF 
Ci=Sxy/Sxs 
PRINT " " 

PRINT USING "12X. ""Sieder-Tate Coefficient = "",D.4D":Ci 
Ac = 0 . 

PRINT " " 

PRINT USING " 1 2X . ""Leas t-Sguares Line: 

PRINT USING "14X. ""Slope = , MD .5DE , " : Ci 

PRINT USING "14X. ""Intercept = "" , MD . 5DE , " ; Ac 

IF Im= 1 THEN 
BEEP 

PRINT USING " 12X . ""NOTE ; "" ,DD. "" data runs were stored in 

ile$ 
run j it 

ASSIGN ®F i le TO * 



file 



• M« 8ft M 



ASSIGN SFilep TO * 

END 

DEF FNRhow(T) 

Ro =1006. 357 24-T*< . 774489-T*<2 . 262459E-2-T*3 . 03304E-4) > 

RETURN Ro 
FNEND 

DEF FNPrw(T) 

P r = 1 0 ‘ ( 1 ,09976605-T*< .01 3759326- T*( 3 . 968875E-5-3 . 45026E-7*T) ) ) 

RETURN Pr 
FNEND 

DEF FNMu(T> 

Mu-1 .57609473E-3-T*(3.51 1 98576E-5-T*(3 . 58358 16E-7-T* 1 .3655861 15E-9) ) 

RETURN Mu 

FNEND 

DEF FNKw(T) 



Nru 



71 



2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 



Ku=. 5721 83504477+1 .527701 21 209E-3-T 

RETURN Ku 

FNEND 

DEF FNTvsv(Emf) 

V=Emf * 1 0'6 

T-V*< .026 173344 16-V*<9. 2447859E-7-V*S . 0746642E- 1 1 >> 

RETURN T 
FNEND 

DEF FNCpu< T > 

Cpu= (4 .211 20853- T*< 2 . 26326E-3-T*< 4 . 4236 1 E-5 + T*2 . 7 1 428E-7 ) ) >* 1 000 

RETURN Cpu 

FNEND 

DEF FNTanh(X) 

P= EXP(X) 

Q=EXP< -X) 

Tanh=(P+Q)/(P-Q) 

RETURN Tanh 
FNEND 
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1000! FILE NAME: PLOT 
1010 PRINTER IS 705 
1020 BEEP 

1 03 0 INPUT "ENTER MINIMUM PND MAXIMUM X-VALUES" , Xmn . Xmax 
1040 BEEP 

1050 INPUT "ENTER MINIMUM AND MAXIMUM Y-VALUES" , Ym in , Ymax 
1060 BEEP 

1070 INPUT "ENTER STEP SIZE FOR X-AXIS" .Xstep 
1080 BEEP 

1080 INPUT "ENTER STEP SIZE FOR Y-AXIS" , Ystep 
1100 BEEP 

1110 PRINT "IN ;SP1 : IP 2000,2000,8000,7000;" 

1120 PRINT "SC 0 , 1 00 , 0 . 1 00 ; TL 2,0;" 

1130 Sf x= ! 00/(Xmax-Xmin) 

1140 Sf y = 1 00/ ( Ymax - Ym m) 

1150 PRINT "PU 0,0 PD" 

1160 FOR Xa=Xmin TO Xmax STEP Xstep 

1170 X= ( Xa-Xmi n ) *Sf x 

1180 PRINT "PA";X,",0; XT;" 

1190 NEXT Xa 

1200 PRINT "PA 100.0;PU;" 

1210 PRINT "PU PA 0,0 PD" 

1220 FOR Ya=Ymin TO Ymax STEP Ystep 
1230 Y=(Ya-Ymin)*Sfy 
1240 PRINT "PA 0 , " ; Y , "YT" 

1250 NEXT Ya 

1250 PRINT "PA 0,100 TL 0 2" 

1270 FOR Xa=Xmin TO Xmax STEP Xstep 

•280 X= ( Xa-Xm in) *$f x 

1290 PRINT "PA" : X , " . 1 00 ; XT" 

1300 NEXT Xa 

1310 PRINT "PA 100.100 P1J PA 100,0 PD" 

1320 FOR Ya=Ymin TO Ymax STEP Ystep 

1330 Y= ( Ya- Ym i n ) *Sf y 

1340 PRINT "PD PA 100,",Y,"YT" 

1350 NEXT Ya 

1360 PRINT "PA 100.100 PU" 

1370 PRINT "PA 0,-2 SR 1 .5,2" 

1380 FOR Xa=Xmin TO Xmax STEP Xstep 
1390 X=<Xa~Xnin)*Sf x 
1400 PRINT "PA" :X , " .0 : " 

1410 PRINT "CP -2,-1 ;LB";Xa;"" 

1420 NEXT Xa 

1430 PRINT "PU PA 0.0" 

1440 FOR Ya=Ymin TO Ymax STEP Ystep 
1450 Y= ( Ya- Ymin) *Sfy 
1460 PRINT "PA 0 . " : Y . "" 

1470 PRINT "CP -4 , - . 25 ; LB" ; Ya; "" 

1480 NEXT Ya 
1490 BEEP 

1500 INPUT "ENTER X-LABEL" .XlabelS 
1510 BEEP 

1520 INPUT "ENTER Y-LABEL" , YlabelS 

1530 PRINT "SR 1.5.2:PU PA 50,-10 CP" : -LEN< XlabelS) /2 : "0 ; LB" : Xlabe IS 
1540 PRINT "PA -11.50 CP 0 . " ; -LEN( Y labels) /2*5/6: "DI 0 . 1 ; LB" ; Y Labe 1$ 
1550 PRINT "CP 0,0" 

1 560 Repeat : ! 

1570 BEEP 
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1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 



INPUT "ENTER THE NAME OF THE DATA FILE" . D_f i leS 
ASSIGN @Fi le TO D_file$ 

BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER". Md 
BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" . Npai rs 
BEEP 

INPUT "SELECT A SYMBOL FOR THE PLOTTER < 1 =* , 2= + .3=c . 4=o . 5=' > " , Sy 
PRINT "PU DI" 

IF Sy = 1 THEM PRINT "SM*" 

IF Sy =2 THEN PRINT "SM+” 

IF Sy =3 THEN PRINT "SMc" 

IF Sv=4 THEN PRINT "SMo" 

IF Sy =5 THEN PRINT "SM'" 

IF Md>Npairs THEN 
FOR 1=1 TO (Md- 1 ) 

ENTER »File;Xa.Ya 
NEXT I 
END IF 

FOR 1=1 TO Npairs 
ENTER ®F i ie ; Xa. Ya 
X= ( Xa-Xmi n ) *Sf x 
Y= ( Ya- Ym i n) *Sf y 
PRINT "PA" ,X , Y , "" 

NEXT I 
BEEP 

ASSIGN SFile TO * 

INPUT "DO YOU HAVE MORE DATA TO BE PLOTTED ( 1 =YES . 0=MO) 7" ,Go_on 
IF Go_on=1 THEN Repeat 
PRINT "PU SM" 

BEEP 

INPUT "DO YOU LIKE TO DRAW A STRAIGHT LINE?".Go_on 

IF Go_on= 1 THEN 

BEEP 

INPUT "ENTER THE SLOPE". SI 
BEEP 

INPUT "ENTER THE INTERCEPT" , Ac 
FOR Xa=Xmin TO Xmax STEP (Xmax-Xmin) 

Ya=Ac+$ i*Xa 
Y=(Ya-Ymin>*Sfy 
X = (Xa-Xtnin)*Sf x 
IF Y<0 THEN 
Xam=(Ynun-Ac)/Sl 
X = (Xa/n-Xmin)*Sfx 
Y=0 

END IF 

IF Y> 1 00 THEN 
Xam = < Ymax-Ac) /S 1 
X=<Xan-Xm in) *Sf x 
Y = 1 00 
END IF 

PRINT "PA" , X , Y , "PD" 

NEXT Xa 
END IF 

PRINT "PU SPO" 

END 
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10 ! FILE NAME; G_L0SS 

20 ! LAST UPDATED: 23 MAY 33 

100 PRINTER IS 701 

110 CLEAR 709 

120 BEEP 

130 INPUT "ENTER MONTH. DATE AND TIME (MM :DD :HH :MM:SS>" .Dates 
140 OUTPUT 709 : "TD" : Dates 

150 OUTPUT 709 ; "TD" 

ISO ENTER 709: Dates 

170 PRINT USING " 1 OX , ""Month . date and tine: "",1 4A" ; DateS 
180 BEEP 

190 INPUT "ENTER INPUT MODE ( 1 =3054A ,2-FILE)" , In 
200 IF Im= 1 THEN 
210 3EEP 

220 INPUT "GIVE A NAME FOR DATA FILE" , D_f i leS 
230 CREATE BDAT D_fileS,15 
240 ELSE 
250 BEEP 

260 INPUT "GIVE THE NAME OF THE EXISTING DATA FILE" , D_f i ie$ 

270 BEEP 

280 INPUT "ENTER THE NUMBER OF RUNS STORED". Nrun 
290 END IF 

300 ASSIGN 3F i 1 e TO D_fileS 
310 BEEP 

320 INPUT "GIVE A NAME FOR THE OUTPUT FILE" ,0_f ileS 

330 CREATE BDAT 0_fileS.15 

340 ASSIGN sFileo TO 0_fileS 

350 K=0 

360 BEEP 

370 INPUT "OK TO GO AHEAD (HIT ENTER)". Ok 

380 IF In= 1 THEN 

390 OUTPUT 709 : "AR AF19 AL21" 

400 OUTPUT 722 :"F1 R1 T1 Z1 FL1" 

410 FOR 1=0 TO 2 
420 OUTPUT 709: "AS SA" 

430 ENTER 722:Enf(I> 

440 IF I>0 THEM Emf ( D=ABS(Emf ( I) >*1 0"6 
450 NEXT I 

460 Tsteam=FNTvsv((Enf( 1 )+Emf(2>)*.5) 

470 OUTPUT 709 : "AR AF23 AL24" 

480 OUTPUT 722: "FI R1 T1 Z1 FLt" 

43U t UK J=U TO 1 

500 OUTPUT 709: "AS SA" 

510 ENTER 722 :Emf ( J+3) 

520 Enf ( J+3 ) =ABS(Emf ( J+3> ) * 1 0"6 

530 NEXT J 

540 OUTPUT §File;Emf(*> 

550 ELSE 

560 ENTER 9File:Emf<*> 

570 END IF 

580 Tsteam=FNTvsv( (Emf(1 > +Emf (2) )* .5) 

590 T room = FNTvsv ((Emf(3) + Emf(4)>*.5> 

600 Pouer=(Enf (0>*100> “2/5 .76*1 .E-3 
610 Dt=Tstean-Troom 

620 PRINT USING " 1 OX , 3D . 3D . 3 ( DDD . DD , 2X ) " ; Power , Tsteam . T room , D t 

630 K=K + 1 

640 IF Irn= 1 THEN 

650 BEEP 
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660 INPUT "ARE YOU TAKING MORE DATA ( ! = YES , 0=N0> . Go on 
670 IF Go_on=1 THEN 380 
680 ELSE 

690 IF K < N r an THEN 380 
700 END IF 
710 END 

720 DEF FNTvsv(V) 

730 T=V»( . 0261 733441 6-V*(9 . 2447859E-7-V*6 .0746642E- 11)) 

740 RETURN T 
750 FNEND 
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Month, date and time: 06:01:14:55:20 

The following analysis was performed for data in file DDtO 



Position number 
Wall temperature (Deg C) 
Average wall temperature 
Inlet water temperature 
Outlet uater temperature 
Log-mean- temp difference 
Water velocity 
Sieder-Tate parameter 
Nusselt number 



12 3 4 5 6 

66.92 65.43 62.00 57.65 65.31 53.03 
61 .56 (Deg C) 

23.94 (Deg C> 

26.32 (Deg C) 

36.42 (Deg C) 

1 .37 (m/S) 

5342.8 

153.72 



Position number 
Wall temperature (Deg C) 
Average wall temperature 
Inlet water temperature 
Outlet water temperature 
Log-mean- temp difference 
Water velocity 
Sieder-Tate parameter 
Nusselt number 



1 2 3 4 5 6 

68.20 66.10 61.57 55.55 65.90 53.93 
61.42 ( Deg C) 

25.36 (Deg C) 

27.39 (Deg C) 

35.03 (Deg C) 

1 .92 (m/S) 

7032.7 

199.76 



Pos i t ion number 
Wall temperature (Deg C) 
Average wall temperature 
Inlet water temperature 
Outlet water temperature 
Log- mean- temp d l f f erence 
Water velocity 
Sieder-Tate parameter 
Nusselt number 



12 3 4 5 6 

68.10 65.35 60.61 54.03 65.49 57.68 
60.43 (Deg C) 

26.77 (Deg C) 

23.59 (Deg C) 

32.79 (Deg C) 

2.47 (m/S) 

8634.8 

247.76 



Position number 
Wail temperature (Deg C) 
Average wall temperature 
Inlet water temperature. 
Outlet water temperature 
Log -mean- temp di f f erence 
Water velocity 
Sieder-fate parameter 
Nusselt number 



:. 1 2 3 4 5 6 

: 67.70 65.11 59.70 52.58 64.40 56.49 

* 59.34 (Deg C) 

= 23.23 (Deg C) 

= 29.94 (Deg C) 

- 30.22 (Deg C) 

= 3.01 (m/S) 

= 10131 .6 

- 302.10 



Position number 
Wall temperature (Deg C> 
Average wall temperature 
Inlet water temperature 
Outlet water temperature 
Log- mean- temp d l f f erence 
Water velocity 
S leder-Tate parameter 
Nusselt number 



12 3 4 5 6 

66.69 64.08 58.48 51.22 63.06 55.10 
58.04 (Deg C) 

29.73 (Deg C) 

31 . 19 (Deg C) 

27.58 (Deg C) 

3.5b < m/S) 

11676.1 

347.19 



Position number 
Wall temperature (Deg C) 
Average wall temperature 
Inlet water temperature 
Outlet water temperature 



1 2 3 

67.18 64.46 53.95 
58.53 (Deg C) 

31 .31 (Deg C> 
33.23 (Deg C) 



4 5 6 

51.63 63.76 55.68 
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Log-mean- temp difference 
Water velocity 
Sieder-Tate parameter 
Nusselt number 



25.93 (Deg C) 
3.84 <m/S> 

1 2518. 8 
371.31 



Position number 
Wall temperature (Deg C> 
Average wall temperature 
Inlet uater temperature 
Outlet water temperature 
Log-mean- temp difference 
Water velocity 
S i eder -Tate parameter 
Nusselt number 



1 2 3 4 5 S 

70.32 G7.3S S2.08 55.23 66.93 59.06 
61 .34 (Deg C> 

33.57 (Deg C) 

35.09 (Deg C) 

27.51 (Deg C) 

3.29 (m/S) 

11223.3 

331.04 



Position number 
Wall temperature (Deg C> 
Average wall temperature 
Inlet *uater temperature 
Outlet uater temperature 
Log -mean- temp difference 
Water velocity 
Sieder-Tate parameter 
Nusselt number 



1 2 3 4 5 6 

74.00 71.28 65.83 58.91 70.63 62.82 
65.70 (Deg C> 

35.46 (Deg C> 

37.16 (Deg C) 

29.38 (Deg C> 

2,74 (m/S) 

9850.0 

284.11 



Position number 
Wall temperature (Deg C) 
Average wall temperature 
Inlet water temperature 
Outlet water temperature 
Log -mean- temp d i f f erence 
Water velocity 
$ l eder -Tate parameter 
Nusselt number 



12 3 4 5 6 

77.43 74.65 69.49 62.86 74.08 66.46 
69.42 (Deg C) 

37.27 (Deg C) 

39, 19 (Deg C) 

31 .18 (Deg C> 

2.19 (m/S) 

3354.3 

240.23 



Position number 
Wall temperature (Deg C) 
Average wall temperature 
Inlet water temperature 
Outlet water temperature 
Log- mean- temp difference 
Water velocity 
S l eder- Tate parameter 
Nusselt number 



1 2 3 4 5 6 

30.70 77.96 73.09 67.18 77.57 70.37 

73.22 (Deg C) 

38.92 (Deg C) 

41 .23 (Deg C) 

33.12 (Deg C) 

1 .64 (m/S) 

6722.9 

200.30 



Sieder-Tate Coefficient = .0294 



Least-Squares 

Slope 

Intercept 



Line: 

= 2 . 93533E-02 

= 0 . 00000E+00 
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Month, date and tune : 05 : 30 : 1 5 : 27 : 1 5 

The following analysis was performed for data in file DR1 



Measured pressure 

Pressure corresponding to measured steam temp 
Noncondensable gas concentration 
Noncondensab le gas concentration 



9 . 578E+04 (Pa) 
9.01 0E+Q4 (Pa) 
5.8 (Mol 7.) 
9.0 (Mass 7.) 



Mater velocity 
Reynolds number 
Inside heat-tran coef 
Qveral 1 -heat- trans coef 
Heat flux 

Outside heat-tran coef 
Nusselt coefficient 
Nus coef (from heat flux) 

Mater velocity 
Reynolds number 
Inside heat-tran coef 
Overal 1 -neat-trans coef 
Heat flux 

Outside heat-tran coef 
Nusselt coefficient 
Nus coef (from heat flux) 

Mater velocity 
Reynolds number 
Inside heat-tran coef 
Overai I-heat-trans coef 
Heat flux 

Outside heat-tran coef 
Nusselt coefficient 
Nus coef (from heat flux) 

Mater velocity 
Reynolds number 
insiae heat-tran coef 
Overal 1-heat-trans coef 
Heat flux 

Outside heat-tran coef 

Nusselt coefficient 

Nus coef (from heat flux) 

Mater velocity 
Reynolds number 
Inside heat-tran coef 
Overal 1-heat-trans coef 
Heat flux 

Outside heat-tran coef 

Nusselt coefficient 

Nus coef (from heat flux) 



= 0.32 (n/S) 

= 13139.2 

= 5356.2 (M/ra'2-K) 

= 3393.6 <W/m'2-K) 

= 2 . 22899E+05 (M/m'2) 
= 19990.5 (W/m'2-K) 

- n304.5 (W/m'2-K) 

= 10797.9 ( W/(tT2-K) 

- 1 .54 (m/S) 

- 26418.9 

= 9074.3 < W/m'2-K) 

= 4109.5 (H/m'2-K) 

= 2 . 63764E + 05 (M/m'2) 
= 11186.9 (W/m~2-K> 

= 10320.1 (H/m'2-K) 

= 1 0005.0 (W/m'2-K) 

= 2.47 (m/S) 

= 40022.4 

= 12466.1 (W/m'2-K) 

= 4725.9 (W/m'2-K) 

= 3 . 0635 1 E+05 (M/m '2) 
= 10550.9 (M/m'2-K) 

- 9822.2 (H/m'2-K) 

= 9434.3 (H/m'2-K) 

= 3.29 (m/S) 

= 53980.1 

= ! 5667 . 3 < W/m'2-K) 

= 5331 .5 (M/m'2-K) 

= 3 . 43200E+05 (M/m'2) 
= 10926.6 (W/m'2-K) 

= 95H.1 (M/m'2-K) 

= 9075.4 (W/m'2-K) 

= 3.34 (m/S) 

= 63586.5 

« 17728.9 (M/m'2-K) 

= 5652.2 (M/m'2-K) 

= 3.61 307E+05 (M/m'2) 
= 11060.4 (M/m'2-K) 

= 9359.5 (M/m'2-K) 

= 8389.1 (M/m'2-K) 



NOTE: 5 X-Y pairs were stored in plot data file PDR1 
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